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^ Abstract: The possibility of testing trilinear neutral gauge boson 

■ couplings in radiative neutral current scattering at HERA is an- 

alyzed using a Monte Carlo program that includes the Standard 
\q ', Model at tree level and the anomalous vertices. Acceptance and 

Q\ \ isolation cuts are applied as well as optimized cuts to enhance the 

signal from new physics. The bounds on Z77 couplings that can be 
achieved are not so stringent as present bounds, even for high lumi- 
nosities, but probe a different kinematical region almost unsensitive 
c-| ! to form factors. 



"*-! ■ 1 Introduction 



The precision data collected to date have confirmed the Standard Model to be 
a good description of physics below the electroweak scale [|X]] . Despite of its 
great success, there are many reasons to believe that some kind of new physics 
must exist. On the other hand, the non-abelian structure of the gauge boson 
self-couplings is still poorly tested and one of the most sensitive probes for new 
physics is provided by the trilinear gauge boson couplings (TGC) 0. 

Many studies have been devoted to the WWj and WW Z couplings. At 
hadron colliders and e + e~ colliders, the present bounds (Tevatron 0]) and 
prospects (LHC, LEP2 and NLC [Q, |J) are mostly based on diboson production 
(WW, W'j and WZ). In ep collisions, HERA could provide further information 
analyzing single W production (ep eWX 0) and radiative charged current 
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scattering (ep — > Z/7X 0). There is also some literature on WW^f couplings in 
VT-pair production at future very high energy photon colliders (bremsstrahlung 
photons in peripheral heavy ion collisions and Compton backscattered laser 
beams ||). 

Only recently, attention has been paid to the Z^Z, Z77 and ZZZ cou- 
plings. There is a detailed analysis of Z7V couplings (V = 7, Z) for hadron 
colliders in ||. CDF and D0 |TT| have obtained bounds on the Z'-fZ 
and Z77 anomalous couplings, while L3 has studied only the first ones ||12| . 
Studies on the sensitivities to these vertices in future e + e" colliders, LEP2 (4| 
and NLC jO, have been performed during the last years. Some proposals 



have been made to probe these neutral boson gauge couplings at future photon 
colliders in e7 — > Ze |T3| . 

In this work we study the prospects for measuring the TGC in the process 
ep — > cyX. In particular, we will concentrate on the Z77 couplings, which can 
be more stringently bounded than the Z^Z ones for this process. 

In Section 2, we present the TGC. The next section deals with the different 
contributions to the process ep — > e'yX and the cuts and methods we have 
employed in our analysis. Section 4 contains our results for the Standard 
Model total cross section and distributions and the estimates of the sensitivity 
of these quantities to the presence of anomalous couplings. Finally, in the last 
section we present our conclusions. 



2 Phenomenological parametrization of the neu- 
tral TGC 

A convenient way to study deviations from the standard model predictions 
consists of considering the most general lagrangian compatible with Lorentz 
invariance, the electromagnetic U(l) gauge symmetry, and other possible gauge 
symmetries. 

For the trilinear Z^V couplings (V = 7, Z) the most general vertex func- 
tion invariant under Lorentz and electromagnetic gauge transformations can be 
described in terms of four independent dimensionless form factors [13], denoted 
by hY, i=l,2,3,4: 



v 



+hle^q 2p + ^ £ ^p p q 2 J. (1) 
z 

Terms proportional to gf and q§ are omitted as long as the scalar com- 
ponents of all three vector bosons can be neglected (whenever they couple to 
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almost massless fermions) or they are zero (on-shell condition for Z or U(l) 
gauge boson character of the photon). The overall factor, f(V), is p 2 — q\ for 
Z'-fZ or p 2 for Z77 and is a result of Bose symmetry and electromagnetic gauge 
invariance. These latter constraints reduce the familiar seven form factors of 
the most general WWV vertex to only these four for the Z'-fV vertex. There 
still remains a global factor that can be fixed, without loss of generality, to 
Qz-yZ = Qz-y-y = e. Combinations of h\{h\) and h\{h\) correspond to electric 
(magnetic) dipole and magnetic (electric) quadrupole transition moments in 
the static limit. 

All the terms are C-odd. The terms proportional to h\ and h\ are CP-odd 
while the other two are CP-even. All the form factors are zero at tree level 
in the Standard Model. At the one-loop level, only the CP-conserving h\ and 
h\ are nonzero [|16j but too small (0(a/ir)) to lead to any observable effect 
at any present or planned experiment. However, larger effects might appear in 
theories or models beyond the Standard Model, for instance when the gauge 
bosons are composite objects HT7H . 

This is a purely phenomenological, model independent parametrization. 
Tree-level unitarity restricts the Z'-fV to the Standard Model values at asym- 
potically high energies [18|]. This implies that the couplings hj have to be 
described by form factors hf(q 2 , q%,p 2 ) which vanish when qf, q\ or p 2 become 
large. In hadron colliders, large values of p 2 = s come into play and the energy 
dependence has to be taken into account, including unknown dumping factors 
@. A scale dependence appears as an additional parameter (the scale of new 
physics, A). Alternatively, one could introduce a set of operators invariant 
under SU(2)xU(l) involving the gauge bosons and/or additional would-be- 
Goldstone bosons and the physical Higgs. Depending on the new physics dy- 
namics, operators with dimension d could be generated at the scale A, with a 
strength which is generally suppressed by factors like (M w / A) d ~ 4 or (y / s/A) d ~ 4 
||19|| . It can be shown that h\ and h% receive contributions from operators of 
dimension > 6 and and h\ from operators of dimension > 8. Unlike hadron 
colliders, in ep —>■ cyX at HERA energies, we can ignore the dependence of 
the form factors on the scale. On the other hand, the anomalous couplings are 
tested in a different kinematical region, which makes their study in this process 
complementary to the ones performed at hadron and lepton colliders. 



3 The process ep — > e^/X 

The process under study is ep — > ejX, which is described in the parton model 
by the radiative neutral current electron-quark and electron-antiquark scatter- 
ing, 

H (-) 

e q -»• e q 7. (2) 
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There are eight Feynman diagrams contributing to this process in the Stan- 
dard Model and three additional ones if one includes anomalous vertices: one 
extra diagram for the Z^Z vertex and two for the Z77 vertex (Fig. 




Figure 1: Feynman diagrams for the process e q — > e 57. 

Diagrams with 7 exchanged in the t-channel are dominant. Nevertheless, 
we consider the whole set of diagrams in the calculation. On the other side, 
u-channel fermion exchange poles appear, in the limit of massless quarks and 
electrons (diagrams (c) and (d)). Since the anomalous diagrams (e) do not 
present such infrared or collinear singularities, it seems appropriate to avoid 
almost on-shell photons exchanged and fermion poles by cutting the transverse 
momenta of the final fermions (electron and jet) to enhance the signal from 
anomalous vertices. Due to the suppression factor coming from Z propagator, 
the anomalous diagrams are more sensitive to Z77 than to Z^Z vertices. In 
the following we will focus our attention on the former. 

The basic variables of the parton level process are five. A suitable choice is: 
(energy of the final photon), cos 7 , cos 6 q i (cosines of the polar angles of the 
photon and the scattered quark defined with respect to the proton direction), 
(the angle between the transverse momenta of the photon and the scattered 
quark in a plane perpendicular to the beam), and a trivial azimuthal angle 
that is integrated out (unpolarized beams). All the variables are referred to 
the laboratory frame. One needs an extra variable, the Bjorken-x, to connect 
the partonic process with the ep process. The phase space integration over 
these six variables is carried out by VEGAS [f2D| and has been cross-checked with 



the RAMBO subroutine [21 



We adopt two kinds of event cuts to constrain conveniently the phase space: 
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• Acceptance and isolation cuts. The former are to exclude phase space 
regions which are not accessible to the detector, because of angular or 
efficiency limitations:^ 

8° < e , 6 ' 7 , 0j et < 172°; 
E e , £ 7 , p% > 10 GeV. (3) 

The latter keep the final photon well separated from both the final elec- 
tron and the jet: 

cos(7, e) < 0.9; 

R > 1.5, (4) 

where R = V Ar] 2 + <f)' 2 is the separation between the photon and the 
jet in the rapidity-azimuthal plane, and (7, e) is the angle between the 
photon and the scattered electron. 

• Cuts for intrinsic background suppression. They consist of strengthening 
some of the previous cuts or adding new ones to enhance the signal of 
the anomalous diagrams against the Standard Model background. 

We have developed a Monte Carlo program for the simulation of the process 
ep — ► e'jX where X is the remnant of the proton plus one jet formed by the 
scattered quark of the subprocess (0). It includes the Standard Model helic- 



itity amplitudes computed using the HELAS subroutines Pz| . We added new 
code to account for the anomalous diagrams. The squares of these anomalous 
amplitudes have been cross-checked with their analytical expressions computed 
using FORM |[23|| . For the parton distribution functions, we employ both the set 1 
of Duke-Owens' parametrizations [25 and the modified MRS (A) parametriza- 



tions |pq| , with the scale chosen to be the hadronic momentum transfer. 

As inputs, we use the beam energies E e = 30 GeV and E p = 820 GeV, the 
Z mass Mz = 91.187 GeV, the weak angle sin^ = 0.2315 [2j] and the fine 



structure constant a = 1/128. A more correct choice would be the running 
fine structure constant with Q 2 as the argument. However, as we are inter- 
ested in large Q 2 events, the value a(M|) is accurate enough for our purposes. 
We consider only the first and second generations of quarks, assumed to be 
massless. 

We start by applying the cuts (0) and (|) and examining the contribution 
to a set of observables of the Standard Model and the anomalous diagrams, 
separately. Next, we select one observable such that, when a cut on it is 
performed, only Standard Model events are mostly eliminated. The procedure 
is repeated with this new cut built in. After several runs, adding new cuts, the 
ratio standard/anomalous cross sections is reduced and hence the sensitivity 
to anomalous couplings is improved. 

2 The threshold for the transverse momentum of the scattered quark ensures that its 
kinematics can be described in terms of a jet. 
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4 Results 



4.1 Observables 

The total cross section of ep — *> e^X can be written as 

O = 0"SM + Ti ' h l + H a i ' ( h D 2 + °12 ' ^1 ^2 + °"34 ' ( 5 ) 

The forthcoming results are obtained using the MRS '95 parametrization of 
the parton densities[] ||26|| . The linear terms of the P- violating couplings and 
hj are negligible, as they mostly arise from the interference of standard model 
diagrams with photon exchange (P-even) and anomalous P-odd diagrams (r 3 ~ 
r 4 ~ 0). Moreover, anomalous diagrams with different P do not interfere either. 
On the other hand, the quadratic terms proportional to (hj) 2 and {hi) 2 have 
identical expressions, and the same for h\ and h\ {p\ = 03, 01 = U4). Only the 
linear terms make their bounds different. The interference terms a 12 and a u 
are also identical. 

We have analyzed the distributions of more than twenty observables in 
the laboratory frame, including the energies, transverse momenta and angular 
distributions of the jet, the photon and the final electron, as well as their 
spatial, polar and azimuthal separations. Also the bjorken-x, the leptonic and 
hadronic momenta transfer and other fractional energies are considered. 

The process of intrinsic background suppression is illustrated by comparing 
Figures ||] and |||. For simplicity, only the most interesting variables are shown: 
the energy £"(7) and transverse momentum Pt{i) of the photon; the angles 
between the photon and the scattered electron (7, e), the photon and the jet 
(7, j), and the scattered electron and the jet (e, j); and the leptonic momentum 
transfer Q 2 (e). In Fig. |2], these variables are plotted with only acceptance and 
isolation cuts implemented. All of them share the property of disposing of a 
range where any anomalous effect is negligible, whereas the contribution to the 
total SM cross section is large. The set of cuts listed below were added to reach 
eventually the distributions of Fig. |3|: 

• The main contribution to the Standard Model cross section comes from 
soft photons with very low transverse momentum. The following cuts 
suppress a 97% of these events, without hardly affecting the anomalous 
diagrams which, conversely, enfavour high energy photons: 

P 7 > 30 GeV 

pi > 20 GeV (6) 
3 The values change ~ 10% when using the (old) Duke-Owens' structure functions. 
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• Another remarkable feature of anomalous diagrams is the very different 
typical momentum transfers. Let's concentrate on the leptonic momen- 
tum transfer, Q 2 = —(p' e —p e ) 2 - The phase space enhances high Q 2 , while 
the photon propagator of the Standard Model diagrams prefer low values 
(above the threshold for electron detectability, Q 2 > 5.8 GeV 2 , with our 
required minimum energy and angle). On the contrary, the anomalous 
diagrams have always a Z propagator which introduces a suppression 
factor of the order of Q\jM\ and makes irrelevant the Q 2 e dependence, 
which is only determined by the phase space. As a consequence, the 
following cut looks appropriate, 

Ql > 1000 GeV 2 (7) 

It is important to notice at this point why usual form factors for the anoma- 
lous couplings can be neglected at HERA. For our process, these form factors 
should be proportional to 1/(1 + Q 2 /A 2 ) n . With the scale of new physics 
A = 500 GeV to 1 TeV, these factors can be taken to be one. This is not the 
case in lepton or hadron high energy colliders where the diboson production in 
the s-channel needs dumping factors 1/(1 + s/A 2 ) n . 

The total cross section for the Standard Model with acceptance and isola- 
tion cuts is o"sm = 21.38 pb and is reduced to 0.37 pb when all the cuts are 
applied, while the quadratic contributions only change from o\ = 2 x 10~ 3 pb, 
a 2 = 1.12 x 10~ 3 pb to a x = 1.58 x lO" 3 pb, o 2 = 1.05 x 10~ 3 pb. The 
linear terms are of importance and change from T\ = 1.18 x 10 -2 pb, r 2 = 
1.27 x 10~ 3 pb to n = 7.13 x 10~ 3 pb, r 2 = 1.26 x 10" 3 pb. Finally, the 
interference term <j\2 — 1-87 x 10~ 3 pb changes to o~ 12 = 1.71 x 10 -3 pb. 

The typical Standard Model events consist of soft and low-p^ photons 
mostly backwards, tending to go in the same direction of the scattered electrons 
(part of them are emitted by the hadronic current in the forward direction), 
close to the required angular separation (~ 30°). The low-p^ jet goes opposite 
to both the photon and the scattered electron, also in the transverse plane. 
On the contrary, the anomalous events have not so soft and high-p^ photons, 
concentrated in the forward region as it the case for the scattered electron and 
the jet. 

4.2 Sensitivity to anomalous couplings 

In order to estimate the sensitivity to anomalous couplings, we consider the x 2 
function. One can define the x 2 , which is related to the likelihood function £, 

as 

x * = -21n£ = 2L ( ff <* - a" + <x°ln^) * L^^f (8) 

where L = N th /a th = N°/a° is the integrated luminosity and N th (N°) is the 
number of theoretical (observed) events. The last line of (S) is a useful and 



7 



familiar approximation, only valid when I a th -a° I /ff° <C 1. This function is a 
measure of the probability that statistical fluctuations can make undistinguis- 
able the observed and the predicted number of events, that is the Standard 
Model prediction. The well known x 2_ CL curve allows us to determine the 
corresponding confidence level (CL). We establish bounds on the anomalous 
couplings by fixing a certain \ 2 — $ 2 an d allowing for the h] values to vary, 
N° = N°(h]). The parameter S is often referred as the number of standard 
deviations or 'sigmas'. A 95% CL corresponds to almost two sigmas (5 = 1.96). 
When cr ~ cr SM + (/ij) 2 crj (case of the CP-odd terms) and the anomalous con- 
tribution is small enough, the upper limits present some useful, approximate 
scaling properties, with the luminosity, 




til(L') ~ \7 — h](L). (9) 



A brief comment on the interpretation of the results is in order. As the cross 
section grows with h] , in the relevant range of values, the N° upper limits can 
be regarded as the lowest number of measured events that would discard the 
Standard Model, or the largest values of hj that could be bounded if no effect is 
observed, with the given CL. This procedure approaches the method of upper 
limits for Poisson processes when the number of events is large (> 10). 

In Fig. |] the sensitivities for different luminosities are shown. Unfortu- 
nately, HERA cannot compete with Tevatron, whose best bounds, reported by 
the D0 collaboration [Ol, are 



\K[\, |^| < 1.9 (3.1), 

\hl\, \hl\ < 0.5 (0.8). (10) 

For the first value it was assumed that only one anomalous coupling contributes 
('axial limits') and for the second there are two couplings contributing ('corre- 
lated limits'). Our results are summarized in Table [|. 

The origin of so poor results is the fact that, unlike diboson production 
at hadron or e + e~ colliders, the anomalous diagrams of ep — > e'yX have a Z 
propagator decreasing their effect. The process ep — * eZX avoids this problem 
thanks to the absence of these propagators: the Standard Model cross section 
is similar to the anomalous one but, as a drawback, they are of the order of 
femtobarns. 



5 Summary and conclusions 

The radiative neutral current process ep — *> e^X at HERA has been studied. 
Realistic cuts have been applied in order to observe a clean signal consisting 
of detectable and well separated electron, photon and jet. 
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HERA 


10 pb- 1 


100 pb- 1 


250 pb- 1 


1 fb- 1 


K{ 


-19.0 
-26.0 


14.5 
19.5 


-11.5 
-16.0 


7.0 
9.5 


-9.5 
-14.0 


5.5 
7.0 


-8.0 
-11.5 


3.5 
4.5 


hi 


-21.5 
-26.0 


20.0 
30.0 


-12.0 
-13.0 


10.0 
18.0 


- 9.5 
-10.0 


8.0 
15.0 


-7.0 
- 7.5 


6.0 
12.0 


hi 


-17.0 
-22.5 


17.0 
22.5 


-9.0 
-12.0 


9.0 
12.0 


-7.5 
-10.0 


7.5 
10.0 


-5.5 
-7.0 


5.5 
7.0 


K 


-20.5 
-27.5 


20.5 
27.5 


-11.0 
-14.5 


11.0 
14.5 


-8.5 
-12.0 


8.5 
12.0 


-6.0 
-8.5 


6.0 
8.5 



Table 1: Axial and correlated limits for the Z77 anomalous couplings at HERA 
with different integrated luminosities and 95% CL. 

The possibility of testing the trilinear neutral gauge boson couplings in this 
process has also been explored. The Z^Z couplings are very suppressed by 
two Z propagators. Only the Z77 couplings have been considered. A Monte 
Carlo program has been developed to account for such anomalous vertex and 
further cuts have been implemented to improve the sensitivity to this source of 
new physics. Our estimates are based on total cross sections since the expected 
number of events is so small that a distribution analysis is not possible. The 
distributions just helped us to find the optimum cuts. Unfortunately, compet- 
itive bounds on these anomalous couplings cannot be achieved at HERA, even 
with the future luminosity upgrades.0 As a counterpart, a different kinematical 
region is explored, in which the form factors can be neglected. 
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Figure 2: Differential cross sections (pb) for the process ep — > e'-fX at HERA, 
with only acceptance and isolation cuts. The solid line is the Standard Model 
contribution and the dash (dot-dash) line correspond to 10000 times the o\ (a 2 ) 
anomalous contributions. 
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Figure 3: Differential cross sections (pb) for the process ep — > cyX at HERA, 
after intrinsic background suppression. The solid line is the Standard Model 
contribution and the dash (dot-dash) line correspond to 500 times the o~\ (a 2 ) 
anomalous contributions. 
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Figure 4: Limit contours for Z77 couplings at HERA with an integrated lumi- 
nosity of 10, 100, 250, 1000 pb- 1 and a 95% CL. 
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